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We experimentally demonstrate an alternative method for the dynamic generation of atomic spin 
squeezing, building on the interplay between linear coupling and nonlinear phase evolution. Since 
the resulting quantum dynamics can be seen as rotation and shear on the generalized Bloch sphere, 
we call this scheme twist-and-turn (TnT). This is closely connected to an underlying instability in 
the classical limit of this system. The short-time evolution of the quantum state is governed by a 
fast initial spreading of the quantum uncertainty in one direction, accompanied by squeezing in the 
orthogonal axis. We hnd an optimal value of = —7.1(3) dB in a single BEC and scalability of the 
squeezing to more than 10^ particles with = — 2.8(4) dB. 

PACS numbers: 03.75.Gg, 42.50.Lc, 03.75.Mn 


Introduction. The efficient generation of highly en¬ 
tangled states is among the biggest challenges in quan¬ 
tum technologies, as they allow approaching the ulti¬ 
mate quantum limits like the elusive Heisenberg limit in 
metrology [T]. Within the last decade, technological de¬ 
velopment has enabled a number of schemes that reliably 
produce entangled many-particle quantum states, rang¬ 
ing from spin squeezed to Dicke states for neutral atoms 
[2HT4]. 

The Lipkin-Meshkov-Click Hamiltonian, originally devel¬ 
oped in nuclear physics m, captures the dynamics of 
interacting particles in two modes. It allows the gener¬ 
ation of a rich variety of spin squeezed m and highly 
entangled non-Gaussian states HZ]. In ultracold gases, 
this Hamiltonian can be implemented using linear Rabi 
coupling and atomic interactions between two internal 
states. With that system, entangled non-Gaussian states 
have been created [18]. In this work, we detail the dy¬ 
namic generation of spin squeezed states. 

The addition of linear coupling to the standard one-axis¬ 
twisting scheme m leads to an exponentially increasing 
quantum uncertainty, implying a corresponding squeez¬ 
ing along the orthogonal axis on a fast time scale. This 
results from the interplay of the twist due to the inter¬ 
action and the rotation (turn) caused by the linear cou¬ 
pling. Semiclassically, this behavior can be understood 
as the dynamics of a classical phase space volume around 
an unstable fixed point. The corresponding emergence of 
correlations has also been studied in a three-mode sce¬ 
nario [TQ[ 1^ . 

Experimental system. In our experiment, the Lipkin- 
Meshkov-Glick Hamiltonian is realized by employing two 
internal states of a ^^Rb Bose-Einstein condensate [21], 
the \D = \F = l^mp = +1) and |t) = |2, —1) hyperfine 
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states of the electronic ground state. We simultaneously 
prepare up to 30 independent condensates, each contain¬ 
ing N = 200 — 600 atoms, in an optical standing wave 
potential m- The large trapping frequencies ensure that 
external dynamics is frozen out. Additionally, the array 
of condensates yields many independent realizations and 
enables scalability to large atom numbers [22] . 

Linear coupling between the two internal states is 





FIG. 1. (color online). Experimental system and classi¬ 
cal phase space, (a) We simultaneously generate around 
30 BECs, each containing N = 200 — 600 atoms, using 
an optical lattice potential. We employ two internal states 
\l) = |E = l,mF = +1) and |t) = |2,— 1) of ®^Rb. The pop¬ 
ulations of the two states are read out by absorption imaging 
after Stern-Gerlach separation, b) The nonlinearity y induces 
an angular velocity which depends on the population imbal¬ 
ance and twists the quantum uncertainty. The squeezing dy¬ 
namics of this one-axis twisting scheme in a single BEC is 
shown on a generalized Bloch sphere (right panel), c) The 
interplay of the nonlinearity and linear coupling (rotation 
around the x-axis with angular velocity Q) creates a twist- 
and-turn scenario. The resulting dynamics features initial 
squeezing and emergence of non-Gaussianity at later times. 
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achieved using a two-photon microwave 6.8 GHz) and 
radio frequency 6.3 MHz) transition. The phase of 
the coupling can be nonadiabatically adjusted by switch¬ 
ing the phase of the radio frequency radiation, which is 
created using an arbitrary waveform generator. In ad¬ 
dition, the coupling power can be altered by changing 
the amplitude of the radio frequency signal and attenu¬ 
ating the microwave by use of a fixed attenuator on an 
RF switch with two ports. The atomic nonlinearity is 
enhanced in the vicinity of an interspecies Feshbach res¬ 
onance between |t) and |^) at a magnetic bias field of 
9.12 G. After each experimental cycle, state-dependent 
detection is implemented using absorption imaging after 
Stern-Gerlach separation of the two components [23] (see 
Fig#). The imaging is performed at low magnetic fields 
1 G) after a ramp-down of the bias field in 300 ms. In 
order to inhibit spin-relaxation loss of the |t) state during 
the ramp-down, we transfer its population to |1, — 1) via 
a microwave tt pulse. 

Theoretical description. In a quantum mechanical de¬ 
scription, the system of N indistinguishable two-level 
bosons can be treated as a pseudospin (J = N/2) and 
displayed on a generalized Bloch sphere. Using the re¬ 
spective creation and annihilation operators of the two 
modes, the ^ component of the pseudospin = |(d|d^ — 

d|d^) = {N^ — Ni)/2 is defined by the population differ¬ 
ence between the two levels, which can be directly de¬ 
tected in the experiment. The orthogonal components 
(coherences) are given by = |(a|d^ + 

Jy = ^(d|d^ — d|d^), fulfilling the angular momentum 

commutation relation [Jj^Jk] = '^^jkiJh 

In this pseudospin picture, our experimental system can 

be described by the Hamiltonian 

n = xJ^,-^Jx + sX, ( 1 ) 

which is a special case of the Lipkin-Meshkov-Glick 
Hamiltonian m- The quantum dynamics of this system 
is governed by the relative strength of the three param¬ 
eters: the nonlinearity y arising from the interparticle 
interaction, coupling strength O given by the microwave 
and radio frequency radiation, and the detuning 6 result¬ 
ing from the mismatch of the coupling and the atomic 
frequency, the AG Zeeman shift as well as the particle 
number dependent mean field shifts. 

The resulting dynamics can be understood from the fact 
that angular momentum operators are the generators of 
rotations. Thus, {xJz)Jz leads to a rotation around the 
z direction whose angular velocity x^z depends on the 
population difference Jz. This can be interpreted as a 
twist (Fig.[^). The one-axis twisting scheme exploits 
this nonlinear twist exclusively. 

The second and third term in the Hamiltonian describe 
linear rotations around the x- and the 2 :-axis, respec¬ 
tively. In the regime A = \Nx/^\ > 1, the relative 
sign between linear coupling Q and nonlinearity y can 
be chosen such that the linear rotation leads to a speed¬ 
up of the shear of the quantum state. This is achieved 




FIG. 2. (color online). Squeezing analysis, a) The gener¬ 
ated states are tomographically analyzed by rotation around 
the x-axis and detection of the population difference Nf — N^. 
The uncertainty of a squeezed state is characterized from the 
fluctuations of repeated experiments using three parameters: 
The extension along the elongated direction (long axis, ^max)? 
the size of the minimal fluctuations (short axis, ^min)? ^^-^d 
the optimal tomography rotation angle amin- b) An exem¬ 
plary tomography result after 15 ms of TnT shows a strong 
modulation of the observed variances and suppression below 
the classical limit (dotted line). 


when the linear coupling transfers the enlarged spread in 
phase (Jx) into an increased spread in particle number 
difference (J^^), which implies faster twisting. 

Classical phase space picture. Further insight can be 
gained by analyzing the corresponding classical descrip¬ 
tion, which is valid in the limit N ^ 00 . The correspond¬ 
ing classical Hamiltonian is 


^class 




z‘^ COS (p + 



( 2 ) 


with the imbalance z = {N^ — N^)/N and the phase 
(p = arctan((J^)/(Ja,)). In the case of sole twisting, the 
phase space portrait features two stable fixed points at 
the north and south pole of the generalized Bloch sphere. 
Exemplary classical trajectories are visualized by solid 
lines in Fig.[^. The addition of linear coupling leads in 
the case of dominating interaction, i.e. A > 1, to two ad¬ 
ditional fixed points on the equator of the Bloch sphere, 
one of which is stable and one unstable imilllS]. The 
TnT scenario exploits this unstable fixed point and is ex¬ 
perimentally realized for A 1.5 and (5 ~ 0 (Fig.[^). 
The instability leads to rapid spreading of the quan¬ 
tum state along the separatrix which divides the classical 
phase space into regions of macroscopically different tem¬ 
poral behavior. Spin squeezing is generated during the 
early dynamics m- At later times, the bending around 
the stable fixed points leads to the appearance of non- 
Gaussianity, and squeezing vanishes imiiH]. 

In contrast, the quantum states created by the one-axis¬ 
twisting scheme HHISIIIII remain Gaussian on much 
longer timescales. 

State preparation. In the experiment, after initial 
preparation of all atoms in 1, a coherent superposi¬ 
tion between |t) and ||) is produced by applying a 7r/2 
pulse of the linear two-photon coupling. Subsequently, 
microwave and radio frequency are attenuated to reach 
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FIG. 3. Spin squeezing dynamics, (a) Spin squeezing for 
different evolution times obtained from one-axis twisting (gray 
squares) and TnT squeezing (black squares), which vanishes 
as the state becomes non-Gaussian. The minimal obtained 
value for the TnT scenario is — 7.1(3) dB after 15 ms evolu¬ 
tion time. The experimental results are in agreement with a 
numerical Monte Garlo wave function analysis for the exper¬ 
imental parameters, which includes the effects of loss (solid 
lines). Additionally including known sources of noise during 
readout yields better quantitative agreement (dashed lines). 
The dotted line depicts the classical limit, (b) In contrast 
to the one-axis twisting scheme (gray), the optimal tomogra¬ 
phy angle for the TnT scheme (black) increases with time, (c) 
The fluctuations along the long axis increase exponentially for 
the first 25 ms in the TnT scheme (black squares, dash-dotted 
line: exponential fit), indicating the underlying classical in¬ 
stability. Error bars are statistical 1 s.d. confidence intervals. 


the regime of A « 1.5 and the phase of the Rabi coupling 
is adjusted by 37r/2, which changes the rotation axis from 
the y to the negative x-axis. Additional phase shifts due 
to the microwave attenuator are compensated by shifting 
the phase of the radio frequency accordingly, and the fre¬ 
quency is adjusted, taking into account the change in AC 
Zeeman shift caused by the altered power of the coupling 
radiation. 

The influence of technical detuning fluctuations, caused 
by variations of the magnetic bias field of 45 pG over 
several days, is reduced by applying a spin-echo pulse 
(tt rotation around the x direction) at half the evolu¬ 
tion time. This also reduces the sensitivity to coupling 
phase errors. By omitting the linear coupling during the 
evolution time, the same sequence is used for a direct 
comparison to the one-axis twisting scenario. 

State analysis. To investigate the states generated by 
the twist-and-turn scheme, we perform a tomographic 
readout. This is achieved by rotating the state around 
the x-direction for various angles and analyzing the fluc¬ 
tuations of the particle number difference for repeated 
measurements (Fig.|^). The fluctuations are quantified 
by the number squeezing parameter = Var(A^^ — 
Ni)/ydiiclass normalized to the binomial variance 
of the corresponding coherent spin state Varciass(^) = 
4p(l — p)N with p = {N^)/N. One representative exper¬ 


imental result is shown in Fig.[^, from which we extract 
the maximal observed variance as well as the minimal 
fluctuations = min(^^) and the corresponding rota¬ 
tion angle (amin- 

Experimental results. The generation of squeezing with 
the TnT scheme is quantified in Fig.[^ by the resulting 
spin squeezing parameter = ^^-^/(cos taking into 
account the reduced mean spin length N{cos (p)/2 due to 
the extension of the state along the long axis. We infer 
{cosp) by applying a rotation with an angle 7r/2 — amin 
and detecting the distribution of particle number differ¬ 
ences. Since the population imbalance after the rotation 
is z = sin p, we can directly access the expectation value 
(cos p)^. 

We experimentally find strong initial spin squeezing, 
reaching a minimal value of = —7.1(3) dB after 15 ms 
evolution time (Fig.|^). The precisely characterized pho¬ 
ton shot noise of the absorption images (standard devi¬ 
ation (Jdet ~ 4 atoms for each component) has been sub¬ 
tracted for all given values. Small fringe noise contribu¬ 
tions remain [23]. Our results confirm the generation of 
entanglement during the early evolution. 

After this minimal value, squeezing is quickly lost. This 
does not imply the loss of entanglement, as spin squeez¬ 
ing only captures the variance properties of the state and 
thus does not fully characterize non-Gaussian states. In 
this regime, entanglement can be shown by extraction of 
the Fisher information m- 

The solid lines in Fig. represent the results of a Monte 
Carlo wave function (MCWF) simulation discussed in de¬ 
tail in the Appendix, and qualitatively agrees with the ex¬ 
perimental data. This simulation includes particle losses 
and the resulting change of the parameters x(A^) and 
S{N) which depend on the atomic density, as well as the 
technical detuning fluctuations. The dashed line includes 
additional known sources of noise during the detection 
process, such as the residual losses due to background 
collisions during the ramp-down of the magnetic field, 
which corresponds to a loss of ^ 8 atoms. 

We also compare this novel scheme with the well- 
established one-axis twisting scheme. While the cor¬ 
responding experimentally obtained value of the best 
squeezing (gray squares in Fig.|^) is comparable, the 
transition to non-Gaussian states happens at much larger 
time scales, which are longer than the largest investigated 
evolution time. 

The experimentally extracted angle of minimal fluctua¬ 
tions Qfmin and the size of the long axis (see Fig. §• and c) 
are much less susceptible to loss than the minimal squeez¬ 
ing, yielding very good agreement with our MCWF sim¬ 
ulations. The underlying classical instability in the TnT 
scheme leads to an exponential growth of quantum me¬ 
chanical uncertainty, which is directly observed by ana¬ 
lyzing the increase of fluctuations of the long axis shown 
in Fig. This is in contrast to one-axis twisting, where 
the deviation from the initial exponential growth occurs 
much earlier m, and confirms that the exponentially 
fast initial squeezing is prolonged in the TnT scheme. 
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FIG. 4. (color online). Scalability to large atom num¬ 
bers. The best obtained number squeezing after 15 ms of 
TnT evolution can be scaled to large particle number by sum¬ 
ming up the atom numbers of adjacent sites (upper left panel), 
yielding = ~3.1(4) dB for 10 200 particles (gray circles). A 
differential analysis of the fluctuations between two different 
parts of the array (lower left panel) is more robust against 
technical fluctuations and improves this value to — 4.5(4) dB. 
The remaining decrease of the observed fluctuation suppres¬ 
sion is caused by the atom number inhomogeneity of the lat¬ 
tice, which leads to different optimal rotation angles. As a 
reference, the scaling for a coherent spin state is indicated as 
black diamonds and is in good agreement with the classical 
limit (dotted line). Error bars are 1 s.d. confidence intervals 
from a resampling analysis [29]. 


nical fluctuations, as these are suppressed for the dif¬ 
ference 5z = 2;ieft — bright [H] . The corresponding rel¬ 
ative squeezing parameter = Var((5^)/Var((5^)ciass 
quantifies the noise suppression relative to the classical 
limit Var((52;)class- With this analysis, we find = 
—4.5(4) dB for the full sample. This is on a comparable 
level with the value of = —5.3(5) dB that can be ob¬ 
tained using one-axis twisting [22] . 

The remaining decrease of squeezing can be attributed 
to the atom number dependence of the TnT squeezing 
scheme (Fig.[^, for which the degree of single-site squeez¬ 
ing improves as the particle number increases. This re¬ 
sults from the particle number dependence of the detun¬ 
ing since the intra-species scattering lengths of the two 
components are slightly different. The numerical solu¬ 
tion of the stationary Gross-Pitaevskii equation for our 
experimental situation reveals to good approximation the 
scalings S oc ^/N and x oc Both effects are in¬ 

cluded in the MCWF simulation (solid lines). Since we 
fix the angle of rotation to a = 52°, the strong depen¬ 
dence amin on atom number (see Fig.[^) mainly lim¬ 
its the scalability using mesoscopic samples with varying 
atom number. Note that this atom number dependence 
of the final state is stronger compared to the one-axis¬ 
twisting scheme, which leads to the slightly lower value 
for the relative squeezing of the full ensemble. 

For the sum of several independent condensates, the two¬ 
mode approximation is not valid, and the mean spin 
length has to be extracted from the visibility V of Ram¬ 
sey fringes. As we observe V = 94.2% for the full en- 


Upscaling of TnT to large atom numbers. The obtain¬ 
able phase precision in an interferometric sequence is lim¬ 
ited by Ac/) = fs/VN for a spin squeezed state with N 
particles and spin squeezing parameter ^s- Therefore, it 
is one of the main challenges to generate squeezing also 
for large atom numbers. It has been shown that the use 
of many mesoscopic condensates allows the upscaling of 
squeezing to large numbers [22] by adding up the popu¬ 
lations of the individual condensates. 

To analyze the number squeezing parameter for different 
ensemble sizes, we sum the atom numbers of different 
lattice sites Vttot = Sj Nfi and V^tot = J2i (upper 
left panel of Fig.|^ and calculate in analogy to the 
evaluation for the single sites. For an evolution time of 
15 ms and the optimal angle a = 52°, this analysis yields 
a noise suppression = —3.1(4) dB even for the full en¬ 
semble of 10 200 particles (gray circles in Fig.[^. As the 
classical reference, the corresponding scaling for a coher¬ 
ent spin state is in agreement with the shot-noise limit 
for all atom numbers (diamonds). 

The loss of squeezing for large atom numbers is caused 
both by atom number inhomogeneities over the array and 
technical noise sources, which are dominated by fluctua¬ 
tions of the magnetic bias field. The influence of the lat¬ 
ter can be minimized by employing a differential analysis, 
in which the array is divided in two parts and the rela¬ 
tive fluctuations of the two population imbalances zieu 
and 2(right are analyzed. This is robust against the tech¬ 




Number of atoms 



FIG. 5. Atom number dependence in the single BECs. 

Due to the criticality of the TnT scheme and the atom number 
dependence of x and (5, the characteristics of final states after 
a fixed evolution time change with atom number. This is 
shown for an evolution time of 15 ms and three key parameters 
of the spin squeezed states: a) The spin squeezing parameter 
^s, b) the optimal rotation angle Omin and c) the extension 
of the long axis. The results are reproduced by a MGWF 
simulation including the experimental parameter dependences 
on atom number (solid line). The dashed line indicates the 
results including additional noise during the readout sequence, 
the dotted lines are the corresponding classical limits. Error 
bars are statistical 1 s.d. confidence intervals. 
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semble, the whole resource can be directly exploited for 
quantum enhanced measurements, either in a DC or a 
gradiometric scheme. The corresponding spin squeezing 
parameters [28] are = — 2.8dB for the direct analysis, 
and ^gdiff = — 4.0 dB for the differential case, showing the 
applicability for quantum-enhanced metrology. 

Conclusion. We have shown that the twist-and-turn 
scheme can efficiently generate spin squeezing on short 
experimental time scales and is thus favorable for squeez¬ 
ing in lossy environments. These spin squeezed states 
can be directly employed in quantum-enhanced measure¬ 
ment schemes, which have already been demonstrated 
for clocks and magnetometry [22j [30H32] . Extending this 
process to ID systems, the criticality of the TnT scheme 
is the underlying mechanism of the miscible-immiscible 
quantum phase transition at zero temperature [33] . Such 
a system is ideally suited for studies of spatial quantum 
correlations and scaling behavior in pattern formation 
dynamics. 
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APPENDIX: MONTE CARLO WAVE FUNCTION 
SIMULATION 

For the quantitative description of the experimental 
results (Fig. 3 and 5), we perform numerical simula¬ 
tions using the Monte Carlo wave function method 
[34H36] . This description includes the effects of atomic 
losses as well as the atom number dependencies of 
the parameters, which change due to varying atomic 
densities. For our calculations, we find good agreement 
using X = 27 r X 1.43Hz/v/]V, ^ = 27 r x 19Hz and 
S = — 27 r X 0.63(\/]V — \/550) Hz. These parameters are 
consistent with independent measurements using plasma, 
TT as well as Rabi oscillations, and the determination 
of collisional shifts from Ramsey sequences. Two-body 
spin-relaxation loss from 11) and three-body Feshbach 
losses are included, with loss rates calibrated by inde¬ 
pendent measurements. The simulation also includes the 
spin-echo pulse at half the evolution time and detuning 
fluctuations with a standard deviation as = 27r x 0.45 Hz 
caused by variations of the bias magnetic field at 9.12 G. 
These fluctuations are independently determined by re¬ 


peated Ramsey measurements. All pulses are simulated 
with a Rabi frequency of flpuise = 27r x 340 Hz and in 
the presence of the atomic nonlinearity, which leads to 
an effective shortening of the spin-echo rotation around 
the X axis. 

For each data point, we numerically calculate 8000 
trajectories using a 4th order Runge-Kutta method 
and evaluate the observables by calculating the mean 
of the expectation values for the different trajectories. 
In each time step of a single trajectory, a random 
number determines if either the wave function is evolved 
according to an effective Hamiltonian incorporating the 
ideal description Eq. 1 with the addition of decay terms, 
or a loss event is implemented |34f[36] . This is done by 
properly cutting and renormalizing the evolved state 
vector and adjusting S{N) and x(^) accordingly. 

In Fig.|^ the results of the MCWF simulation are 
compared with the ideal theory. We assume an initial 
atom number of N = 500 atoms and the time-averaged 
parameter Nx = 27r x 30 Hz, J = 0 Hz and U = 27r x 19 Hz 
for the twist-and-turn scheme, and U = 0 Hz for one-axis 
twisting. On a qualitative level, we find good agreement 
between the ideal quantum evolution and the results 
of the MCWF simulations shown in Fig. 3 of the main 
text. For the chosen parameters, the value for the 
optimal spin squeezing obtained from the ideal evolution 
is = —I3.0dB after an evolution time of 18 ms, while 
for the MCWF description we find optimal squeezing of 
= —lO.IdB after 16 ms. 




FIG. 6. Comparison of ideal theory and MCWF re¬ 
sults. The results of the MCWF simulations (thin black line 
for the twist-and-turn scheme and thin gray line for one-axis 
twisting) qualitatively agree with the ideal quantum evolution 
(thick solid lines) which does not include losses and the corre¬ 
sponding parameter changes for (a) spin squeezing parameter 
^s, (b) optimal rotation angle Omin, and (c) the fluctuations 
along the long axis of the state. In both ideal theory and 
MCWF simulation, the spin squeezing rate of the twist-and- 
turn scheme exceeds the corresponding rate of one-axis twist¬ 
ing at intermediate evolution times, and squeezing is lost in 
the later evolution. 
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